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We discuss the capabilities of spherical antennae as single multifrequency detectors of gravitational 
waves. A first order theory allows us to evaluate the coupled spectrum and the resonators readouts 
when the first and the second quadrupole bare sphere frequencies are simultaneously selected for 
layout tuning. We stress the existence of non-tuning influences in the second mode coupling causing 
draggs in the frequency splittings. These URF effects are relevant to a correct physical description 
of resonant spheres, still more if operating as multifrequency appliances like our PHCA proposal. 

PACS numbers: 04.80.Nn, 95.55.Ym 



A spherical antenna is particularly well adapted to 
sense metric tidal gravitational wave excitations as a con- 
sequence of the perfect matching between the GW ra- 
diation patterns and the sphere's vibration eigenmodes. 
When suitably monitored, it can act as an individual 
multimode detector generating potential knowledge on 
the incident direction of the signal, its quadrupole or 
monopole amplitudes, or its polarizations |]l],D. 

Spherical geometry also offers another substantial abil- 
ity which optimizes the capabilities of resonant bars cur- 
rently in operation By comparison, on equality 
of building material and mass, the estimated sphere's 
absorption cross section for the fundamental mode is 
over a factor of 4 better than that of a Weber's cylin- 
der in its most favorable orientation. In addition, the 
second sphere's quadrupole still shows a rather large 
cross section, only about half its value at the lowest fre- 
quency and still about 15 times bigger than that associ- 
ated to the cylinder for the fundamental mode, which is 
the only showing a significant cross section in bars 
Hence, spheres show good sensitivity not only at the first 
but also at the second quadrupole harmonic, and there- 
fore can efficiently operate at two frequencies as single 
multifrequency detectors. 

The challenge becomes how to take advantage of this 
potentiality. A preliminary task before multifrequency 
implementation design is to acquire a clear understand- 
ing of the detector's performance when working at one 
resonant frequency -see for instance 1^ or ^ . 

The theories are proposed for practical multimode de- 
vices consisting of a perfect elastic sphere of radius R, 
mass M, density g and Lame coefficients A and fj,, mon- 
itored by a suitable readout system, commonly a set 
of J identical resonators modelled as linear springs of 
mass m = riM,ri << 1, and endowed at certain po- 
sitions rio = Xa/R on the sphere's surface. Each res- 
onator amplifies the elastic radial deformation Ua{t) of 
the sphere at its attachment point caused by an exter- 
nal force and supplies an essential increase in the cou- 
pling by measuring the springs' deformations qait) rela- 



tive to the sphere's undeformed surface. Resonant tun- 
ing requires the sensors to be typically built to possess a 
natural frequency f2 equal to one of the eigenfrequencies 
of the free sphere's spectrum, say watl- This is always 
taken to be a monopole or a quadrupole frequency -ujno 
or UJN2 respectively-, since only monopole or quadrupole 
sphere's spheroidal modes can be excited by an incoming 
wave 0. 

In the universally assumed ideal approximation, the 
tuning frequency wjvl is thought to be an isolated res- 
onance frequency (IRF): no other frequency Uni of the 
free sphere's spectrum interferes in the resonance. For 
instance, this assumption is correct when the resonators 
couple to Wi2, in practice the most interesting coupling 
for unifrequency quadrupole radiation sensing |Q] . 

We follow the perturbative approach in |6j naturally 
opened to refined analysis, and recall that the responses 
of the coupled device to an incident GW can be inferred 
from the system 



dt^ 



HV^Ua - (A + A*) V(V-Ua) fr, 



qait)+nlqait) = -Ua{t)+Ca,Gwit) fl = 1, . . . , J, (1) 

where fres{xa,t) and fcwixajt) are respectively due to 
the resonators' back action on the sphere and to the inci- 
dent GW which also causes a tidal acceleration Ca,Gw{t) 
on resonator a. 

After implementation of Green function formalism, the 
equations can be rewriten in the s-Laplace domain as 



J 

E 
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Mabqb{s) 
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Ca,Gwis) 
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with Ua,Gw{s) the bare sphere's radial responses to the 
signal at the resonators' locations. 

In the IRF circumstance, matrix Mab is of the form 



Sab + rj- 



P. 



L.ab 



(3) 
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The geometric properties of a particular resonator ar- 
rangement are displayed by the J x J symmetric real 
matrix Pl associated to the L-mode selected for tuning. 
It basically has as element ab the Legendre polynomial 
of order L, a sum of products of spherical harmonics: 



L.ab 



\ANLiR)\' 



m—L 

J2 Yl^jina)YLM{nh) (4) 

m— — L 



with An l {R) radial components in the spheroidal normal 
modes of the free sphere . It is highly remarkable that 
all the information determining the distinctive readout of 
a given configuration is just concentrated the eigenvalues 
and eigenvectors of Pl- 

First order resolution for the spectrum of coupled- 
mode eigenfrequencies from det Mab — shows that the 
attachment of resonators causes the tuning frequency 
lonl to split into J symmetric pairs around the original 
value, 



1 ± ^a^r- 



oiv) 



1, 



,J, (5) 



whereas resonators amplitudes amount to be linear com- 
binations of the GW amplitudes g'^^"^\s): 



m=2 



X E >'2™(nf,)g(^™)(s) +0(0) a - 1, . . . , J. (6) 

m=-2 J 

Fl(±6) = ±[aNLANLiR)] (-!)■' (Cc)"' 

with unl non zero overlap coefficients only if L = or 

L = 2a. 



Layout 


IRF Eigei 
around u)i2 


frequencies 

around 0^22 


PHCA 


1 

ujo,± = 1-^12(1 ± 0.58r;2 ) 
1 

uji^± = aji2(l ± 0.88r;2 ) 
i^2,± = '^'12(1 ± 1.07r?i) 


iJo,± = 1^^22(1 ± 0.03r?^) 
1 

ciJi,± = tiJ22(l ± 0.05r;2 ) 
uj2,± = '^'22(1 ± 0.06r;i) 


TIGA 


L0± = aji2(l ± I.OO775) 


1 

iO± = tJ22(l ± 0.05772 ) 



TABLE I. IRF frequencies around 0^12 and LU22 for the 
PHCA and TIGA arrangements. 

However, this is not always the case that the inter- 
esting frequency for tuning is an isolated resonant fre- 
quency. For example, the second quadrupolc LU22 is in 



fact a suitable sphere's frequency for a second resonator 
set to be tuned to in order to exploit the antenna as 
a multifrequency device. A careful examination of the 
spectrum of a typical planned full scale aluminium sphere 
-rjs « 1/40000, R = 1.5m- around UJ22 shows that UJ14 is 
merely within a distance of order 77 2 : 



■^14 



^22 



1 + if 77s 



(7) 



where the dimensionless parameter K takes the negative 
value K = —2.21 for rj = rjs. 

The expectation is that this nearness alters in some 
way the above IRF results. The effects of such unisolation 
-URF effects must be accurately considered to reach 
a faithful description of the detector's real behaviour. 

Let us analyse the general case when cjjvl is the single 
selected frequency for tuning and lun'L' is in its neigh- 
bourhood. Again, we start from expressions (1) and (2) 
which remains unchanged, although a new contribution 
appears in (3): 
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We can still go further in drawing generic conclusions 
valid for any resonator distribution whenever it allows 
Pl and Pl' to commute: [Pl,Pl'] = 0. Then, we name 
ti("' the orthogonal basis which simmultaneously diago- 
nalises the two matrices, with eigenvalues £,a l l' 
respectively, so that in this basis Mab is also diagonal. 
Then, the URF resonances around lvnl can be written 
in J triplets 



0(77) 

1 = 1,. 



(9) 



where for each resonator index a, the upper label {T} 
represents the group {u, c, d} which refers to the three 
different solutions of the cubic equation 



X-Ua + 



Kea.2 
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(10) 



with parameters which are univocally determined by fix- 
ing the layout, lonl and uj^il'- Inspection of orders 
of magnitude in (10) for ujnl = ^22 and lun'L' — ^14 
leads to the conclusion that the triplets present a general 
pattern with independence of the resonator distribution: 
one of the three frequencies will always be located very 
close to the original tuning frequency 0^22 -to assess how 
much near this frequency actually is one needs to restrict 
to particular cases and numerical evaluations-. The re- 
maining pair forms a non-symmetric doublet around it. 
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in good agreement with the idea that the presence of 1^14 
causes a perturbation of the ideal IRF situation. 

Then, the URF effect results in a dragging effect break- 
ing the symmetry of the IRF doublets which approach the 
disturbing frequency tuu, and moreover induces the ap- 
pearance of a third central component near the resonant 

0J22- 

Near actually means really near, at least for two known 
proposals: PHCA [|jlO) and TIG A |Jll[l. It can be seen 
from the numbers in Table II that in each a-group the 
central is itself of order 772 or smaller, so that these 
central URF resonances actually differ from W22 in terms 
of order 77. Reproduction of the amplitudes evaluation 
process demostrates that the contribution of such modes 
are not at leading order 77^2 but at higher terms. There- 
fore, they are referred to as be weakly coupled. 



Layout 


Strongly Coupled 


Weakly Coupled 




1 

too,u = i^22{l + 0.06r?2 ) 
t^o.d = ^^22(1 - 1.16772 ) 


U22{1 - 3.3 ■ 10-^77*) 


PHCA 


toi,u = 1^22(1 + lAlrji) 
ioi,d = ^^22(1 - 2.51772) 


0^22(1 - 1.5 ■ 10"''7;i) 




UJ2.U = ^1^22(1 + 1.077?^) 

1 

i^2,d = ^^22(1 - 2.I8772 ) 


L02,c — 

U22{1 - 3.5 ■ 10-^77^) 


TIGA 


LJl_5,„ = LJ22 (1 + 1.2277^) 
LJl_5,d = tJ22(l - 2. 32772) 


1.^1-5, c = 

0^22(1 - 2.5 • 10-''7;^) 




tJ6,S = t^22(l — 1.105775) 


t^e.-D = <^22 



TABLE II. SCD+WCS URF triplets around LJ22. Calcula- 
tion have been performed for the proposals PHCA and TIGA, 
and for 773, a theoretical value for a full scale future spherical 
detector. Subindex u (up) labels the values which are above 
the resonance frequency tJ22, whereas d (down) labels those 
underneath, and c (central) those practically at U22- 

The result is that both PHCA and TIGA URF triplets 
are composed of a a weakly coupled singlet plus a strongly 
coupled doublet (triplet named of the SCD-I-WCS type) 
of span comparable to that of the IRF pairs. The only 
exception is the sixth mode in TIGA which corresponds 
to a triplet composed of a strongly coupled singlet plus 
a weakly coupled doublet (SCS-I-WCD type). Also, the 
degeneracy pattern of the IRF triplets is maintained: for 
PHCA, three different triplets, two of them doubly de- 
generated, and for TIGA one triplet five-fold degenerate 
plus one single weakly coupled triplet. 

As said, only strongly coupled frequencies (SC) con- 



tribute to the amplitudes: 

'?a(s)=??-^EE|^^i'(f^' 
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with frequency weights 
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(12) 



The high degree of symmetry showed by both PHCA 
and TIGA explains why these configurations fulfil the 
property [Pl,Pl'] = for L = 2 and L' = 4, so that 
P2 and P4 simultaneously diagonalise. After algebraic 
calculus, the eigenvectors of the common basis are finally 
found to be arrangements of spherical harmonics of order 
two, precisely the eigenvectors of P2 when diagonalised 
separately in the IRF situation: 



(13) 



By inserting these functions in (12), we arrive to simpli- 
fied expressions for the readouts of the PHCA layout: 



qa.,PHCA{s) = rj 2 

+ ^24(C/f^) 



SC I 



y2o(n.).9(^°)(.) 
(s'+^o,sc) 



+ F2i{use) 



(s^+^?,sc) 

l22(na)ff(22)(s) + Y2-2{na)g'^^-^\s)\ 



(s2 + ) 



(14) 



and the readouts of the TIGA layout: 

m=2 



J2 Y2m{na)g^'"'H-^) 



qa,TIGA{s) = ?7 ^ XI -^24(C/ff 



SC \ m=-2 



SC 



(^'+^?~5,SC) 



(15) 



The weight functions i^24(t^c)^'^ in (14) and (15) can be 
settled from (12). It is here important to recall that 
the overlap coefficient is zero, which is in accordance 
with the idea that only monopole or quadrupole sphere's 
modes can be excited by an impinging GW. One has to 
take also into account that just up and down frequencies 
in the triplets correspond to SC values, and that the the 
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remaining give place to weakly coupled resonances 
which do not contribute at leading order. Hence, 



1 



22 



u.d TjUjd ' 



(16) 



In Figure 1 we display numerical simulations for the 
IRF and URF qualitative responses caused by the ac- 
tion of a GW burst travelling down the symmetry axis of 
the PHCA pentagonal layout. In both cases the practi- 
cal outputs become beats, although differing in modula- 
tion and modulated frequencies. The IRF output exactly 
cancels at the nodes because the wheights affecting the 
contributions to the amplitudes of lo+ and lo- equal each 
other whereas the URF readout shows a certain thickness 
at those same points caused by the difference between 
the wheights associated to lOu and cod- The occurrence 
of different wheights is in direct relation to the dragging 
induced by the URF effect which breaks the symmetry of 
the IRF-doublets. Nevertheless, the similarity of the pat- 
terns again corroborates the weakness of the new third 
central frequency in the URF triplets. 
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FIG. 1. PHCA beat outputs for GW burst with h+ = 
propagating along the pentagonal symmetry axis at t=0. 



Hence, the analysis of the URF effect is essential for 
a complete study of any spherical multifrequency GW 
detector, but we are specially interested in it with respect 
to our PHCA proposal. 

It is actually idealized as a multifrequency antenna 
with two sets of supplementary pentagonal layouts, one 
tuned to the first quadrupole harmonic Wi2 and another 
coupled to the second quadrupole harmonic UJ22, which 
in fact has ujn very close to it. The pentagonal hex- 
acontahedron polyhedric shape with more than enough 
number of faces for 10 resonators matching in non paral- 
lel pentagonal positions (and even an eleventh resonator 
position for monopole sensing) guarantees technical fea- 
sibility, also supported by the absence of cross interac- 
tions between the outputs of each set at order 77^ 2 . Our 




FIG. 2. PHCA Proposal. Marked faces in the pentagonal 
hexacontahedron indicate resonator positions: a square for 
transducers tuned to LJ12, a triangle for that tuned to ^1^22, 
and a star for the monopole sensor. 



model ensures a correct description of the multifrequen- 
cial PHCA coupled spectrum. We find two main groups 
of frequencies, the first composed of pairs symmetrically 
distributed around UJ12 and the second arising as a non- 
symmetrical splitting of LU22 in triplets dragged towards 
UJ14. 
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FIG. 3. PHCA multifrequency spectrum. 

The PHCA monitored at their first two quadrupole 
modes can for instance be advantageously used to de- 
tect chirp signals from coalescencing binary systems, and 
even to determine some of their characteristic parameters 
by means of a robust double passage method [|l2|. This 
potentiality, unthinkable for currently operating bars, 
shows that resonant spherical antennae are abreast of 
projected broadband large laser interferometers with re- 
spect to their predicted ability in monitoring gravita- 
tional waves from these sources . 

In any case, the URF effect plays per se an essential 
part in a rigorous theoretical analysis of any single mul- 
tifrequency GW detector. Its features are concisely re- 
ported by our developments without severely complicat- 
ing the evaluations. We conclude by emphasizing that 
the philosophy underlying these algorithms is easily ex- 
tensible to the study of other practical real situations 
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departing from ideality. 
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